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Abstract. Using a fermionic version of the Lax pair formulation, we construct an integrable 
small-polaron model with general open boundary conditions. The Lax pair and the boundary 
supermatrices K± for the model are obtained. This provides a direct proof of the integrability 
of the model. 
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1 Introduction 

Recently, there has been renewed interest in integrable boundary systems |], ||] due 
to their connection to the Kondo problem || and boundary conformal field theory H) 
in low-dimensional quantum many-body systems. Sklyanin jHj proposed a systematic 
approach to construct and solve integrable quantum spin systems with open boundary 
conditions (BC). Central to his method is the introduction of the so-called reflection 
equations (RE), which are the boundary analogues of the Yang-Baxter equations 
(YBE) ||. For integrable fermion systems, in particular models of strongly correlated 
electrons, graded versions of the YBE and the RE fi| have been applied. However, 
solutions to the graded RE have so far been found only for certain models due to 
the mathematical complexity of the graded RE. An alternative method to Sklyanin's 
approach is the Lax pair formulation for integrable quantum systems with open BC 
JO. In this article, we use the Lax pair formulation to construct an integrable model 
of correlated fermions with general BC. The Lax pair and the boundary supermatrices 
K± are given explicitly. 

2 Fermionic version of the Lax pair formulation 

We first recall the Lax pair formulation for integrable fermion models with general 
open BC in one dimension We consider an operator version of the auxiliary linear 
problem 

= Lj(u)9j, j = l,2,...,N, 



5*' 



with boundary equations 



dt 



A/ y (>/.)<!>.,, J = 2, 3, 



= M+(u)$jv+i 



,iV, 



(1) 



dt 



$j = M_(u)$i. 



(2) 



Here, Lj(u), Mj(u), and M±(u) are matrices depending on the spectral parameter it; 
u itself does neither depend on the time t, nor on dynamical variables. Evidently, the 
consistency conditions for Eqs. ([!]) and (||) yield the Lax equations 



-L 3 (u) = Af i+1 («)i 3 -(u) - L^uJM^u), j =2,3,. 



,JV-1, 



(3) 



and 



M + (u)£jv(u) - LjvWMjvW, 
M 2 {u)L 1 {u) - Li(u)M_(u). 



(4) 



Let T(u) — Ljv(u) • • • L2(u)Li(u) be the usual monodromy matrix || of the system 
and if_ (it) , K + (u) the supermatrices for the left and the right boundary, respectively. 
The transfer matrix t(u) is defined as the supertrace on the auxiliary space Vq as 



-(u) = str [if + (u)T(u)tf_(u)T -1 (-iO] 



(5) 



From the Lax equations (|^) and (Q), it follows that the transfer matrix t(u) does not 
depend on time t, provided the boundary matrices satisfy the conditions 



str 



K-(u)M-(-u) 
K + (u)M+(u)f(u) 



M-(u)K^{u) 



= str 



K + (u)T(u)M+(-u) 



(G) 



where T(u) = T(u)K-(u)T 1 (—u). This implies that the system possesses an infinite 
number of independent conserved quantities, and thus is completely integrable. 



3 The small-polaron model with general open BC 

We consider the small-polaron model H , which describes the motion of an additional 
electron in a polar crystal with external field parallel to the transverse direction. The 
Hamiltonian reads 

N N N 

h = -jY,^ a i-i +a \-i a i) +v Y. n i n i-i +w Y. n i 

j=2 j=2 4=1 

+ p + n N + a+ajy + (3 + a N + p^n 1 + a_a| + (3^a 1 , (7) 

where in the boundary terms p±, a± and /3± are Grassmann variables, with p± 
even and a± , {3± odd. Hermiticity of the Hamiltonian requires «j_ = j3± and p± = 



2 



p± The fermionic creation and annihilation operators a- and satisfy the usual 

anticommutation relations, and rij = a\a,y We parametrize the coupling parameters 
J, V, and W as 

J = 1, 

7 = -2cos(2?7), 

W = 2sin(2ri)tan(w) + cos(2r/), (8) 

in terms of rj and w. The monodromy matrix corresponding to the Hamiltonian (Q) 
is given in Ref. Jl(|, and obeys the graded YBE. This ensures the integrability of 
the model with periodic BC. In what follows, we shall use this monodromy matrix 
and construct the boundary if-matrices as described in Sec. ^| in order to show the 
integrability of the Hamiltonian. 

4 Constructing the Lax pair at the boundaries 

The equations of motion —i-^O = [H,0] for the fermionic operators at the left 
boundary are 



1 — Hi 

dt 1 
d 

- l dt^ 



dt 



-Ja\ + Vn 2 a[ + (W + p-)a\ + 
Ja 2 — Vn 2 a 1 — (W + P-)a 1 + a- , 
— J(a\a 1 + a 2 a\) — ct-a\ + (3-a 1 . 



(9) 



At the right boundary, they are given by 



dt 

d 
i — 

dt 

d 
i—' 
dt 



N 



'TV 



'N 



-JOn-1 + Vn N~l^N + ( W +P+) a N + P+i 
Ja>N-l - Vn N~l a N - ( W + P+) a N + a +- 



J( a N a N-l + a N a N-l) - a + a N + P+ ( 



■N- 



(10) 



Using the Lax matrices Lj(u) and Mj(u) for the bulk part as given in Ref. |U 
together with Eq. (|j), we can construct the boundary Lax matrices 



(11) 



After some algebra, we find 
1 



A±(u) 



D±(u) 



S+2S-2 

+s T 2 (is±2 =F CJ 2 sinw) (3±a± 

1 

S+2S-2 

+s±2 (is T 2 ± C+ 2 sinu) j3±a ± 



i sin 2 2ij (p± + cos 27/) (1 — n ± ) + s±2 (iszp2 ± £+ 2 sinu) a±a± 

isin2?7tanLj, (12) 
— i sin 2 2r] (p± — cos 2r/)n :j _ + s T 2 (is±2 T £+ 2 sin u) a±aj_ 



isin2?7tana;, 



(13) 



3 



B±{u) 



C±(u) 



S+2S-2 



— C+ 1 sin u sin 4?? =p sin2?7 (p± sinw ± sin2?7COSM)a_(_ 



sin 2?7 s T2 a± 



(14) 



sin u sin Arj f3±n± — z£+ 1 sin 277 (p± sin u =p sin 2r\ cos it) 



S+2S-2 

±^i si n2?7 S±2 /3 ± 

For convenience, we introduced the abbreviations 
s± n — sin(u ± nrj), 

a_|_ = djy, 

and 



c±„ = cos(w ± n?7), 



COs(li ± bS) 



(15) 



cos it cos u> 



(16) 
(17) 



5 Boundary reflection supermatrices 

We are now ready to compute the boundary A-matrices 



K±(u) 



-"■21 -"22 



(18) 



by using Eq. (^|). This yields 16 equations for the four matrix elements i£y, only six 
of which are independent. After a lengthy calculation, we find 



- [1 n — 

AT, = 



A 10 = 



sin(u - %/}-) (C-S+2 - C+s-2) , 



a_ sin"0_ sinu , 2 



sin 277 



^21 = 



/3_sin^_sinu , 2 2 w 2 2 



-"■li — 

#2+2 = 



i sin 2?7 

sin(u + 2?7 - (C-S+i - £+ sinu) 
£_ s'm(u + 2r\ + tp + ) (£ + s+4 - £_ sin it) 



K^2 = 



Ktx = - 



a + sinijj+s+2 
sin 2 2<q 
/3 + sin-0+s +2 / 2 



(£ 2 s+4 ~ £+ sinu) (C+s+4 - £ 2 sinw) 



i sin 2?7 



(C-S+4 - £+ sinu) (C+s+4 - £ 2 sim 



(19) 
(20) 

(21) 

(22) 

(23) 
(24) 

(25) 
(26) 



where we have parametrized p± = sin 2r] cot i/j± . The Hamiltonian (Q) can be found as 
usual as an invariant of the commuting family of transfer matrices t(u) (|^) by taking 
the derivative at a special value of the spectral parameter u. Namely, 



sin2r7 — r(u) 
du 



u=0 



= 2ffr(0) + str -rK + {u) 
1 du 



u=0 



(27) 
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6 Conclusion 



We have presented the Lax pair formulation for the ID small-polaron model with 
general boundary conditions, providing a direct demonstration for the integrability 
of the model. The Grassmann variables at the boundary terms play a crucial role 
in the construction of the boundary if-supermatrices. They may be interpreted in 
physical terms as sources and sinks for injecting additional particles into the system. 
In particular, we can choose as a special representation at the left end a_ = c_a , 
ft- = c_ Oq, and p_ = c_ , and at the right end a+ = c+a N+1 , /3+ = c^a) N+l , 
and p + — c?) , with complex numbers c± and two additional boundary fermions a 
and a N+1 . We note that the boundaries considered here act as pure back-scatterers 
[[[], p). Thus a combination with the forward-scattering impurities of Refs. [pjj |l2| 
may result in physically relevant, yet completely integrable impurities. Corresponding 
Hamiltonians may be constructed by the methods outlined here |D| . 
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